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Abstract 

A spherically symmetric evolution model of self-gravitating matter with the equation of state p = —{1 + 6)e 

(where 6 = const) is considered. The equations of the model are written in the frame of reference 
r^ I comoving with matter. A criterion for the existence and formation of a horizon is defined. Part of 

^P ■ the Einstein equations is integrated analytically. The initial conditions and the constraints imposed 

on these conditions in the presence of a horizon are determined. For small 6 an analytic solution 
, to spherically symmetric time-dependent Einstein equations is obtained. Conditions are determined 

Q ■ under which the dynamics of matter changes from collapse to expansion. Characteristic times of the 

^^ . evolution of the system are evaluated. It is proved that the accretion of phantom matter (for 5 > 0) 

onto a black hole leads to the decreases of the horizon radius of the black hole (i.e., the black hole is 

"dissolved"). 
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7 1 INTRODUCTION 

Matter that violates the zero energy condition (NEC) is said to be phantom matter |lj. For matter 
that violates this condition, the sum of the energy density e and pressure p, (e + p) is less than zero. 
^ ■ In the present paper, we assume that the energy density e of phantom matter is positive (in the 

^O , co-moving frame of reference). The sum {e + p) is proportional to the energy density in the frame 

^^ I of reference co- moving with a photon p]. Of course, such a frame of reference is senseless; only the 

(^ • asymptotic limit for a frame of reference whose velocity approaches the velocity of light makes sense. 

Phantom matter has a number of exotic properties; this fact gives a motivation for its comprehen- 
sive study. 

First, only in this matter "wormholes" can exist (see, for example, [3j or [3j). 
Second, there is a conjecture [5j that phantom matter possesses a unique property to dissolve black 
holes in it. Until now, this has been proved only for a non-self-consistent solution (when the effect 
p\ ' of the gravitation of this matter can be neglected compared with the gravitation of ordinary matter), 

j^ . The complete dissolution of a black hole in phantom matter would require revision of the definition 

and the description of the term " event horizon" (see Section [3|) . 

Third, the state of the art in observational cosmology allows one to speak of possible prevalence 
of phantom matter in the universe (according to the latest measurements of the acceleration of the 
expanding universe). In addition, a real example of phantom matter is given by the quantum field of 
the Casimir effect (between close conducting plates [6j). Thus, the study of such matter is of practical 
interest [7], [3]. 

Fourth, there is an analytic solution to the Einstein equations for the cosmological equation of state 
e + p = 0, therefore, small deviations from this exact equality are possible toward the equation of state 
corresponding to the case of phantom matter, and it becomes possible to find a first-approximation 
solution in a small parameter. 

The analytical solution thus obtained removes the above-mentioned contradictions and questions. 
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2 EQUATIONS OF MOTION 

Let us choose a frame of reference co-movmg with matter (see, for example [8]). It is convenient to 
choose the metric tensor in the spherically symmetric case aqj: 

ds^ = e'^dt'^ - e^dR^ - e^'idO^ + sin^ Odf^). (1) 

Here, the metric component e'^ = r^ determines the area of the sphere around the center of the system 
(47rr^); the metric components r, u and A are functions of the coordinates R and t. 
The Einstein equations corresponding to metric ([!]) can be expressed ao 

Stte = -e-^ (/x" + f (/i')2 - i/i'A') + e-- (i/iA + ^fif) + e'^ , (2) 

Sttpii = -e-'^ (/i + |(A)2 - i/iz>) + e-^ (i/xV' + \{fi'f) - e'^ , (3) 

Svrpx = e-^ i^^v" + i(i/')' + y + 3(/^')' " 3/^'^' " i^'^' + i/zV) + 
+e- (1 Az> + i/,z> - 1 AA - iA - i(A)2 - i/i - i(/i)2) , 

2/*' + /i^u' - A^' - i/'/t = . (5) 



(4) 



where py = — T^ is longitudinal pressure, p_L = — T^ = —T]^ is transverse pressure, the prime denotes 
the derivative with respect to iZ, the dot denotes the derivative with respect to t, and Tf are the 
components of the energy- momentum tensor. 

If pressure is isotropic, py = p_l = p, then these equations yield two useful relations, which can also 
be obtained directly from the formula Tf.^ = (the Bianchi identities): 

A + 2/1 = -2e/(p + e) , v' = -2p'/{p + e) . (6) 

Let us write out the equation of state of the matter: 

P = P±=P\\=we, w = -(l + 5). (7) 

where tu is a constant parameter that determines this equation. When 5 = 0, Eqs. ^ clearly show 
that —p = £ = const, this implies the well-known solution [9j with a A-term: 

£ = e^ > , v^ = -\^ = ln(l - rg/R - AR^) . (8) 

In this expression, the constant r^ corresponds to the gravitational radius in the presence of a horizon, 
while the constant e^ = 3A/(87r) is the energy density of the cosmological A-term. 

Due to the presence of the A-term, the metric becomes unphysical for large R (due to the external 
horizon). Therefore, here one should define the maximum possible radius Roo'- 

ARl = l. (9) 

In addition, the asymptotics of the metric for large R must be Galilean: 

lyoc^O- (10) 

Below, when speaking of physical infinity, we will assume that the boundaries of the metagalaxy are 
not reached, i.e., R < Roo- 



^We use a system of units in which c = 1 and G = 1. 

^These equations are derived, for example, in ^ (Problem 5 in §100). 
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r < Th 




Figure 1: The straight hne a — b — c — d — eis the world line of an observer in the coordinates of the 
CO- moving frame of reference. The curves represent the graphs of t(/j) for various surfaces r = const. 
The dashed lines are light cones within which light can propagate. At the point b, the light cone 
intersects the surfaces r = const; hence, there is no horizon at this point. At the point d, the light 
cone passes above the curve r = const; hence, it intersects only the surfaces with r < r/^ j); therefore, 
the point d is below the horizon. 



3 CRITERION FOR THE EXISTENCE OF AN APPARENT HORI- 
ZON 

The radius of the event horizon depends on the entire future evolution of matter [3] . This is inconve- 
nient for describing real black holes from both mathematical and practical points of view because the 
definition of the event horizon is not local. Usually, when speaking of a black hole, we assume that 
there exists an apparent horizon and that the event horizon will be formed at the end of evolution (it 
will be clear from the sequel that this does not always correspond to reality). 

Let us give a local definition of the apparent horizon rh [3] . We prove that the apparent horizon in 
a spherically symmetric system is formed at the moment when an incident particle with nonzero rest 
mass reaches the velocity of light with respect to the surfaces r = const at the same point at which 
the radiating particle is located at this moment. 

Now suppose that we are on a particle of matter with coordinate r^o and, from a large radius Roo, 
follow up a particle of matter with coordinate Rp that radiates light while passing through spheres 
of radius rrji ^^y A criterion that a particle has not reached the horizon is the fact that we can still 
see light emitted by this particle; i.e., light intersects surfaces with radii r = const. Hence, a criterion 
that a particle reaches the horizon is the event that the propagating light cannot intersect surfaces 



with radii r > r 



{Rp,t)- 



Let us express this criterion mathematically. 



In Fig. [H the straight vertical line a — b — c — d — e denotes the world line of the particle Rp in 
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the coordinates of the co-moving frame of reference {R, t) from the moment of rest (a) to the center of 
the system (e), where r = 0. The sohd curves passing through the points e, d, c, and b represent the 
hues of constant values of ^(/j j) for r = 0, r < r/^, r = r/j and r > rh, respectively. The dashed lines 
with vertices at these points indicate cones within which light emitted by particle Rp can propagate. 
Therefore, according to the above criterion, the horizon exists at the point where the cone is tangent 
to the line r = const; in Fig. [H this line corresponds to the value r = r^ (it passes through the point 
c). 

The condition of constancy of the radius rm^\ has the form: 



dr = = rdt + r' dR. 
Hence, the tangent of the slope of the curve r = const with respect to the axis R is 



dt 
dR 



r= const 



r 

r 



For the light cone, we have (by definition) ds'^ = 0, which, with regard to ([T|), yields 



dt 
dR 



v^- 



ds=0 



Then, a criterion for the absence of a horizon is given by the condition: 

dt 



dt 
dR 



ds=o dR 



(11) 



(12) 



(13) 



(14) 



r= const 



Substituting (|T2|1 and (|T3[) into this expression and introducing the physical velocity V of the particle 
of matter, we rewrite this criterion as 



V' 



X-u 



/\2 



{r/r'Y <1. 



Introduce an element of proper (physical) time 

dr = exp{i'/2)dt . 
a physical longitudinal distance / and an element of physical length dl: 



R 



kR,t) = J e^'\R,t) dR, dl = exp(A/2) dR . 



(15) 



(16) 



(17) 



Then, the velocity of motion of matter with respect to the surfaces r = const is given by 



V 



dl 



V^~ 



dR 



Ve^-^- 



(18) 



This proves the assertioiu that the apparent horizon is formed at the moment when matter reaches 
the velocity |F| = 1 with respect to the surface r = const. 



■^This assertion is inapplicable at the mouth of a wormhole (if any) because e '*' identically vanishes there. Therefore, 



the existence condition of a horizon near the mouth is determined asymptotically: V 



1 at r — *■ rji. 
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4 GENERAL INTEGRALS OF MOTION 

Integrating the first equation in ^ with regard to the initial conditions, we obtain 

The integration of the second equation in ([6]) yields 

1 + 6 



6 



ln(e/ej + c„(t), (20) 



By an admissible transformation of time Co(t), we can set the function Co(t) equal to zero. Then, 

u = -2l±^ln{e/ej (21) 

Then, condition (jlOp is fulfilled automatically as i? ^ oo. 
Let us rewrite Eq. ([5]) as 

2r' = r'X + ri^' . (22) 

Determine A from this equation and substitute it into Eq. ^. We obtain 

„2 _ 1 ^ f^ ^^'2^_A , ■2-u 



Sirer^ = — r - rr'e'^ + rr^e'" . (23) 

r an V / 

Let us introduce a new quantity that has the meaning of mass: 

R 

M^R^t^^Jiner^'dR, (24) 



Upon integrating Eq. (I23p . we can obtain two equations 

2M = r(l-r'V^ + rV-) , v"" = 1 - ^ ■ (l -^V (25) 

Hence, provided that r' e~^ ^ 0, the horizon radius is 

r=rh 

TK = 2M^RM ■ (26) 

Note that Eq. ^ coincides with Eq. ([2]) if we replace the time derivative by the derivative with 
respect to the coordinate (and vice versa), make the change X ^^ u, and change the signs in the last 
terms of these equations. Then, by analogy, Eq. (|23p can be rewritten as 

— Svrpii r = ~~r y ~ ''"^' e~ + rr e~'^] . (27) 

Integrating with respect to time and taking into account (j25p and the initial conditions, we obtain 

t 

„2 



2M = 2Mo - / 8TTp\\r^rdt , (28) 



Here, 

Mo^M(^,o), V) = 2M„/i?. (29) 

Formulas (j24p and (|25p yield an expression for the initial mass in the presence of a horizon: 

R 
2Mo =ri+ I 8iT£oR^ dR , (30) 

ri 

where r^ is the initial radius of the horizon, which is defined by r^ = rg + Arj* ~ rg. Formula (j28p 
determines a shift of the horizon in matter due to the work of pressure forces on the internal layers of 
matter (when matter is collapsed to the center). In [lOj, formula (j28p was obtained for a particular 
case oi w = const. 
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5 INITIAL CONDITIONS IN THE ABSENCE OF A HORIZON 

Let us specify the initial conditions of the problem for a model in which there is no horizon at the 
initial moment. Using an admissible transformation R = R(j^\, we can specify: 



,=R =^ r'=l. (31) 



at the initial moment. In addition, we can specify the condition that matter is at rest at the initial 
moment: 

fi,=0, K=0, X,=0, z>„=0. (32) 

Then, Eqs. i^5^ yield 

\ = -ln(l-5), (33) 

where 

S = 2MJR, 

this result does not depend on the choice of 6. 

6 INITIAL CONDITIONS IN THE PRESENCE OF A HORIZON 

In this model, it is assumed that a horizon may exist at the initial moment; this assumption also 
concerns a purely vacuum solution with a A-term. Therefore, the initial conditions near the horizon 
deserve special mentioning. 

The initial condition (|3ip is retained. With regard to (jlSp . on the one hand, we have 

yoiR) = ^'^of ■ exp(Ao - i/JI^^ = 1 , 
while, on the other hand, 



0, e^o 



-^ oo . 



Therefore, we obtain indeterminacy of the form • oo on the horizon. Removing this indeterminacy, 
we obtain a physical distribution for velocity V^, at the initial moment: 



\Vo{R)\ 



» (^^ •■'•!■ (34) 

1 {R = rh) . 



i.e., the initial distribution of velocity represents a function with a "punctured point." 

To validate formula ()34p . we modify the initial conditions: we assume that r^, ^ 0. Rewrite the 
second expression in ([25|) as 

e-^o =i + a-S, (35) 

where a = r'^e~'^o is determined by the initial distribution of the velocity of matter. 
Hence we obtain the following expression for the initial distribution of velocity: 

K' = r—^ . (36) 

" I + a — S 



As a ^ 0, formula (|36|) turns into (p^ . 
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7 ANALYTICAL SOLUTION 

Denote 

y = \n{e/e^). 

As (5 ^ 0, we have e ^ e^; therefore, y ^ as (5 — > 0. 

Substituting A and v' from Eqs. ([6]) into Eq. (122^ . we obtain 



(37) 



^, r' 2 



y 

r 



:i+'^)7 



(38) 



Let us pass from the old coordinates (R, t) to new coordinates (r, r^), which are functions of the old 
coordinates. Applying a trivial mathematical transformation, we obtain 



y y 



dy_ fr-h rh' 
dvh 



(39) 



In Section [TOl we will show that, for small 5^ the variation of r^ is also small (and proportional to 5). 
Taking into account that y — > as (^ — > 0, we can see that the right-hand side of ([39|) can be neglected 
as 5 ^ 0. Then, Eq. (f38l) is rewritten as 



^, r' 2 



0. 



(40) 



Now, let us find a solution to the Einstein equations in which the function y (as well as e and v) 
depend explicitly only on the functions r and r^ and do not depend explicitly on t and R. When the 
density of matter is defined so, the function e(r) is determined solely by the initial distribution of the 
density of matter. 

As (5 — > 0, expressions ([2T]) and ([35]) for u^ and Aq must tend to Uj^ and A^, respectively [see dH)], 
and a must tend to zero. Accordingly, we choose the energy density distribution as 



e = e. 



1 - I In (1 + a 



rh' 



z/=-2i±^y«(l+5)ln(l+a-^) . 
In the main approximation in 5 with regard to the initial conditions, we obtain 



(41) 



2r 
r 



r ~ — 



(42) 



Expression ([30j) yields 



2Mo « rg + Ai?3 _ I^Ar-S . I{R/r^ - 1) 



I{x) = x^ 



\nx 1 

3 9 



+ X^ 



\nx 1 

2 4 



+ X [In x — 1] — (x + 1) 



ln(z+l) _ 1 
3 9 



9 ' 



(43) 



lim I{x) = 3;(ln X — 1) — > , 
X— >o 



where r ^ = fg + Kr^ ~ rg . 

Taking into account that e depends explicitly only on r, we obtain the following expression from 
dSlD: 

2M = 2Mo {R) - (1 + 5) / Svrer^ dr = 2M^ {r) - 6 J Svrer^ dr . (44) 
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In the main approximation in 6, formula (I44p is rewritten as 



2M « Tg + Ar" 



6A 



R' 



r^ + 77^^ • I{r/r^ - 1) 



From (fTOl) and (l35ll we find 



(i?/r) 



1 + a-S 



1 - I In (1 + a - rh/r) 
l-U^C^+a-rh/R) 



2/5 



Hence, using (j25p and applying the second remarkable limit to (06]), we obtain 



(R/rY 



1 + a- Th/R 



V' 



1 



1 - 2M/r fl + a- rn/R 



(45) 



(46) 



(47) 



1 + a — S \1 + a — rfi/r J 1 + a — S \1 + a — r^/r 

According to the second expression in (j25p . the horizon radius rh is determined by the formula 
r/i = 2M. Then, from (|45p we find that 

rh-r^^ -5K{R^ - rf) . (48) 

as (5 — > 0. However, the solution obtained does not clearly indicate the direction of evolution (either 
collapsing or expansion). The direction of evolution is determined by the acceleration of matter at 
the initial moment. 



8 ACCELERATION 

Equations ([3|) and (j35p yield the following expression for the second derivative of radius r with respect 
to time t at the initial moment: 

2Re-''or, = (1 + S)8TTeR^ + Rr^v^e'^o -S + Ru',{l + a - S) . (49) 

At the initial moment, the proper acceleration a = dV/dr of the observer co-moving with matter is 
given by (see (fT8|) ): 

Oq =r; •exp(Ao/2-zyJ. (50) 

For small 5, from (1411) and (j43l) we obtain 

u^^{l + 5)\n{l + a-S) , 



S^^ + AR^- |(5A-i • I{R/r, - 1 



RS' w -^ + 2Ai?2 + |5a5 • [I{R/r, - 1) - RI\R/r, - 1) 
Let us rewrite (I49p with regard to ()51|1t 



2Re-''or^ ^6 



LA + aR^ + -Ar3 . /' 



3 

— j_. 
2 1 



+ CK- 



j/R - 2AR^ + fJArf (I' - I/R) 



1 + a-S 

Hence, we can see that the acceleration vanishes (as it must) at 5 = (and a = 0). 
For radii R close to fj^ = rg + Ar^, the terms containing 

/' « ln(i?/r, - l)/r, < 0, 

prove to be dominant in (j52p . Denote 



/? PS r 



1 + exp 



3(5Ar^i?c 



1 + exp 



36Ar^J 



(51) 



(52) 



(53) 



Then, from (j52p we find that, for 5 > 0, matter expands at the initial moment when R > Rcr and 
collapses when (R < Rcr) (in the close vicinity of the horizon). 



'Here, as before, the prime denotes diflFerentiation with respect to R. 
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10 10 
Figure 2: The behavior of a cloud of phantom matter. The cloud, which has the form of a Gaussian 
distribution when at rest, starts do expand with time. The vertical axis represents the energy density 
normalized to unity. 

9 CHARACTERISTIC EVOLUTION TIMES OF THE SYSTEM 



In this section, we will not assume that S is small. Then, neglecting a and using 
the following expression for the characteristic evolution time of the system: 



91), we can obtain 



T 



r/\r\. 



If the initial distribution of matter is characterized by large gradients, one can neglect all the terms 
in (j49p except the terms containing Ri^g' . Then, (j49p and (j50p yield 



w 



r,^j{e,'/e,){l-S){eJe,r/', 



w 



%- ■j(eo7eo)Vl-5. 



(54) 



According to formulas ()54p . there are two alternative versions of the evolution of matter with the 
initial distribution in the form of a spherical layer with a local maximum of energy density: 

1. In the case of phantom matter (5 > 0), we obtain a expanding layer of matter (see Fig. [2]). 

2. When 6 < 0, we obtain a collapsing layer of matter. 

Note that the greater the density gradient, the faster the evolution. Apparently, these results can 
be generalized to a non-spherical model. 

The characteristic evolution times corresponding to expressions (I54p are given by 



^1 ~ ^\{6/w)Re,(eJe,y-'/s\ / |e/(l - S)\ , T2 ^ ^\{6/w)Re,\ /\e,'/{ 



[1-S) 



(55) 



Here, Ti is the evolution time for an infinitely distant observer, and T2 is the evolution time for an 
observer that moves with matter. 
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To conclude this section, we consider the evolution of the domain that is characterized by the 
maximum energy density (a "hump"; see Fig. [2|). In this domain, Rv^^' ~ 0; therefore, formulas (|49p 
and (150p are rewritten in this domain as 

f, ^ exp(z.„) [(1 + 5) ■ Stte.R^ - S] /{2R) , 

a,^[{l + 6)-8TTe,R^-S]/[2R^/T^^]. ^ ' 

For the initial distribution of matter similar to that shown in Fig. [2l we have 

JSTreR^dR 

S = '^ < Sne.R^, 

R 

at the maximum; therefore, the acceleration at this point proves to be positive; i.e., the "hump" 
decomposes and flies away from the center of the system (for any 5 > 0). 

10 DYNAMICS OF THE APPARENT HORIZON 

Since all particles may intersect the horizon only in one direction, a strict correspondence is established 
between the coordinates Rh and th on the horizon. This means that Rh{t) is a continuous and 
monotonically increasing function {Rh > 0). 

Consider formula (j28p on the horizon. Differentiating this formula with respect to time, we obtain 
the following expression for the variation of the horizon radius: 

rh = SiTEgRh ■ —. Snwer'^r . (57) 

Taking into account (fT2]l . we rewrite (|57|l as 

rh = Sne.Rl ■ Rh ■ ll + wr' ■ ^ ■ ^\ (58) 

Consider two ways in which the horizon radius varies. 

10.1 Test Matter 

Consider the dynamics of the horizon for test matter (which negligibly affects the variations of the 
system and the radius of the horizon). For such matter, all functions must depend only on r. Other- 
wise, the metric in a static (fixed) frame of reference turns out to be nonstatic, which implies that it 
is affected by matter (i.e., the assumption that it is a test matter is violated). 

Equation (jlip determines the constancy of r, while the constancy of the metric component u is 
given by 

du = = u'dR + i>dt. (59) 



Solving Eqs. pTj) and (|59|) simultaneously, we obtain 

ly'r = z>r' . (60) 

Equations ([22]) and ([60]) yield 



2r'/r' = A + z> . (61) 

Integrating this equation with regard to the initial conditions, we obtain 



r' = A/exp (A - Aq + 1/ - i/J . (62) 



10 
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Hence, taking into account (jlOh and (I2ip . we find 

r' = ^~. (63) 

Now, formula (j58p for the variation of the horizon radius is rewritten as 

rh = 87re,Rl-Rh-{l + w). (64) 

Thus, for w < —1, we have rh < 0; i.e., in the case of test matter, the horizon radius decreases, and 
the black hole is "dissolved" in phantom matter. 

In formula (j64|) , the fact that matter is test manifests itself in that the quantity SttEq R^ is extremely 
small, St^EqR? <<< 1. This result was predicted and proved for the static frame of reference (a frame 
of reference fixed with respect to distant stars) in [S]. 

10.2 Nontest Matter 

Formula (|64p can be generalized to nontest matter. To this end, it suffices to notice that, when 
the horizon radius is constant, the metric coefficients on the horizon are also time-independent (the 
gravitational field at radius r is determined only by the mass within this radius, i.e., by the mass 
within the horizon radius in the case considered). Therefore, the following relation must hold for 

dl'\rh = 0. 

Using this in expression (j59p on the horizon, we obtain the same result (j64p . Hence, if the accretion 
of matter does not change the horizon radius, then the equality w = —1 must hold. 

du\^^ = dt + v' dR . (65) 

Solving Eqs. (fTTI) and (i65]) simultaneously, we obtain 

(66) 



/• ■ / t dv 
V r = vr — r — 

which is an analog of Eq. ([60|) . Equations ([66]) and ([22]) yield 

2r'/r' = A + z> - — 

ax j,_,^^_j,Q^g^ 

Integration of this equation with the initial conditions yields analog of Eq. ([62 



(67) 



r' = ^ /exp(A - \ + v - v^) ■ exp [v^^ - z^[2] ] , 



Here, the following notation is used: [1] = (r/i,0) and [2] = {R^^th); r^i^ = r^j = r^. An expression for 
the variation of the horizon radius is obtained analogously (I64|) : 

rh = Stte.RI -Rh-il + w^], ^ = ^exp (z.[i] - z^pj) = ( ^ J • (69) 

These are exact expression independent of 5. Thus, the direction in which the horizon is shifted 
depends on the initial conditions of the distribution of matter. The condition under which the horizon 
radius decreases is given by u;^ < — 1. 

For a matter distribution that depends explicitly only on r, we obtain ^ = 1. Then, formulas (I13p . 
([HD, dUD and ([MD yield 

This formula shows that, in the general case, the "dissolution" degree of a black hole for small 6 is 
rather small. 

11 
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11 CONCLUSIONS 



The mam conclusion of the present study is the confirmation of the conjecture on a possible decrease 
of the horizon radius (solution) of a black hole. 

Note that, at infinity, there is no paradox related to this phenomenon. In the absence of matter, 
the metric at infinity persists to be Schwarzschild (with a A-term). Accordingly, the total mass of 
the whole system related to this metric remains constant. The point is that the mass and energy are 
redistributed in space during evolution: an external wave of phantom matter removes the mass and 
energy that correspond to the initial black hole (which is "dissolved" under the action of the internal 
wave of phantom matter). If the black hole is not "dissolved" completely, the total energy of the 
system is composed of the mass of the changed black hole and the energy of the external wave of 
phantom matter. 
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